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1. INTRODUCTION 


The theoretical determination of the electron current (J.) of the metallic elec- 
trode, placed in ionised gas is an important problem of the plasma physics. А 
dependence of the current on the electron energy distribution [F(s)] is specially 
important. The widely known Druyvestein theory [1] is not fully satisfactory for 
resolving the problem. This theory does not explain the fact, observed in the expe- 


riment with the electrostatic probe, that electron current density 1 = F decreases 


when the area surface S of the probe increases, e.g. see Herrmann and Klagge [2]. 
Attempts at an explanation of the experimental fact were often undertaken. The 
consideration of Swift [3] concerning the spherical probe are successful. However, 
the results of Swift theory can be used only for a low pressure with regards to 
accepted assumption. У 

In the present paper, the generalization of the Swift theory is described. In 
these consideration the problem of the calculation of the special integrals [9(Р, z)] 
is present. Values of many different integrals can be determined by integral 6 (P, z). 
For the above reasons, integral is treated as the new special function. The values 
of 0(Р, z) calculated with the computer are presented in the tabele. 
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2. GENERAL ANALYSIS OF THE PROBLEM 


The metallic electrode, placed in the plasma, is the object of an intensive neu- 
tralization of charge carriers: the plasma becomes locally impoverished. When the 
electrode potential (V) in relation to the plasma is negative, the positive charge 
sheat surrounds the electrode. This sheat having the size (a) shields the electri- 
cal field of the electrode. The determination of the d is an open problem, e.g. 
see report [4]. In general, it is accepted that d is equal to the Debeye length 
um EA _ Gor" f 

у V пее(Ть + T;) 
electron and ion temperatures, n is the plasma density, e is electron charge. If the 
mean free path (Л, } of the electrons is longer than the sheat size, the Druyvestein 
theory (modified by Medicus [5]) can be used for the determination of the electron 
current of the electrode. According to the above theory, the group F1(s)A«s of the 
electrons, placed at distance of d from the electrode, forms the current AJ, defined 
as 


, Со is here the perinitivity of the free space, Т.Т; are the 


д9 о alt= £ File) Be, (1) 


5 
Ге : 
where g = 4/ —— and m is here the electron mass. 


m 
The advantage of the Druyvestein theory is, that it gives a simple relationship 


f a eles, 
of Е, (=) with the second derivative —— of the electron current, 


ау? 
42.7. 4. 95 Е, (e) d (2) 
dV? a 


The above formula is widely used in order to determine the electron energy di- 
stribution function Fy(¢) of an undisturbed plasma. Using the above formula one 
accepts that Ё (=) = Fo(s). However, this assumption is often incorrect. This 
assumption is insignificant if one uses the Druyvestein theory in the Swift modi- 
fication. According to Swift, the plasma perturbation around the electrode has a 
more complex structure than that accepted in Druyvestein theory, see Figure 1. 
Besides the region of the strong electrical field , the region of the weak field is 
present at which the plasma density is lower than that of the undisturbed plasma. 
The electron current passes through this second region in a diffusive manner. At 
distance of г from the spherical electrode, the group |F(s) As] of the electrons forms 
the electron current 
d Fle)Ae| 


dr ‘ (3) 


AJ, = 4rr72D 


; 1 | ZEE Р 4 
where D — 3 My X In case of the spherical electrode, a very useful assumption 


can be accepted that Ё(=) — К (=) as r — oc. When we use the above assumption 
then resolving equation (3), we obtain 
AJ. 


4neD(a + d) ` (4), 


Е (=) = Fo(e} == 
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region of strong Неа 


м. 
Ы? д undisturbed plasma 
FIE ) &F,(£) 


Fig. 1. Plasma stucture arround the spherical electrode 


Substituting (4) into formula (1) and integrating over region V + co we obtain 


JV) = gS(P — 1) [= 25 VePolelde (5) 


cn ee 3a? К | ^ 
here Р = ——, р = -————— / 
where 24 j' № la+ d) and a is the electrode radius. On the basis of 
formula (5) a very important relationship of Fo(e) ee 
d^g, F(e ) 2 [^ Fleve 
уз = 95 fet 2P(1- P) bh e- Py?" (6) 


Swift obtained analogous expressions. The presence of the integrals in formulae 
(5, 6) causes that the further solution of the problem is very difficult. Swift con- 


E UN А didi 
sidering the relationship р 


the Maxwellian and Druyvestcinian distributions of the electrons. Swift calculated 


= f{F(e)] determined the correction function for 


: : А а 
suitable integrals in а rough manner, when he assumed that — < 0.6. However, 
- Ld 
this assumption limits the use of the Swift solution in the experiment with the 
electrostatic probe. 


The typical energy distributions [F (e) $ of the plasma electron can be described 
by the formula: 


F(s) = M БЕ)! (7) 


where a, 8, К are the suitable constants dependent on the distribution type. Sub- 
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stituting (7) into (5, 6) we obtain 


JV) = srol n} [ e - V(1- Р)+ 
mv | (8) 
emn Les]. 


я eco f 1 
—2Р(1-Р) n |. 


(8) |а)" 


2 : 
In case of the Maxwellian distribution (k = 1, a = "d В = 1), the integrals 


d? J. gSna pv j 
ave VTS {52-х 
+ 2PV P?Vv? 

(= - PV)? | fe = PV. 


of formula (8, 9) can be simply determined with the exponential-integral function 
oo 


Е: (т) = у : -exp —t- dt see ref. [6]. For example the: J, (V) can be described as 


X(V) = 95" VT( - P) (oo - 7.) m 


-P(- P) [exw - 14 . E; [a E 24 (10) 


3. CASE OF THE DRUYVESTEINIAN ENERGY DISTRIBUTION OF ELECTRONS 


For stable homogenous plasma (a model plasma), kept with the stable electrical 
field, the Druyvesteinian type of the electron energy distribution is characteristic. 
This distribution is the result of the solution of Boltzmann equation, if the inde- 
pendence of the electron collisions cross section on the electron energy is assumed, 
see Rutsher [7]. It is characteristic of the classical Druyvesteinian distribution that 

zre(2) 
- = 0.565. Often the simple Druyve- 


3975 (4) 


steinian distribution is used for which k 2, В = 0.5, a = - = 0.577 аге 


- 08 
| у21(2) 
characteristics Г(2) = [У t?~1 - exp —t dt is the gamma function, see ref. [6]. 

In case of the Druyvesteinian distribution the integrals of the expression (8, 9) 
cannot be determined by special functions known so far. However, they can be de- 


termined by the expression containing the error function erfc (z) — Tos exp —t? dt 
со 


see ref. |6], and the new special function O(P, х) = [ D -exp —t? dt. 
- Рх . 


т 
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Je Е 
уз can be described as 


According to above, the formulae for J, and for 


J.(V) И СЕ ет] + 


» (11) 
-(1- Р), erfc (вх) - P(1- DAY , 
23. = th —4P? + ЗР? + 2P*(1 - P)(8— =). [exp - (87 298 
—4Р?(1 – PPPZ = 287.) erfc (y (12) 


-2P|(1 - P)? - s(1 ~ PP(87. + РУУ) US Bn} 


ES UN E 


The manner of reducing of integral [ (2 Pv ds to integral 0(P, г) is descri- 


bed in Appendix. 


4. CALCULATION OF Diu) 6(P, x) 


Introducing new variable z — F integral (Р, 5 can be described as follows: 


pix ехр-— 


% exp —t? 
қа) = f vr P wc 


This transformation permits to calculate value 9 (Р, z) by function RUMOAA (dou- 
ble precision) which is described in the reports [8, 9]. The use of the above function 
allows to calculate 0(P, z) with relative errors < 1: 10-9. The calculations were 
made with the computer R-22 (я IBM-360). The computed values of 6(P, =) for 
P=0.1,0.2,...0.9 and z = 0.1, 0.2, ...3.0 are presented in Table 1. Moreover, 
Figure 2 presents the dependence of 0(P, z)on P and т. 

In order to determine intermediate values 0;(P,,z;) (for z;, P, absent in Ta- 
ble 1) linear approximation cannot be used, because the errors can be higher by 
several per cent. The shape of diagrams 0(Р, т) = f(z) shows that the depen- 
dence can bé approximated exponential type functions: 6(P,z) = Aexp—b,/z, 
(Р, =) = Aexp—bz, 0(P,z) = Aexp—bz,/z, at intervals 0.05 + 0.3, 0.2 + 0.8. 
0.7 + 3, respectively. When constants A and b are determined on the basis of 
01(P, 21) and 65(P, zz) values, corresponding to zı , 12 = гу + 0.1 arguments, then 
intermediate 9, (Р, z;), corresponding to zè, (zi < Z; < 12) and determined accor- 
ding to the above approximation, is charged with an error smaller than 0.25%. In 
case of the calculation of intermed'ate values 8; (Р., т) (for intermediate P, parame- 


can be 


ters, different from those given in Table 1) approximation 6(P, z) = 
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used. When constants A and b are calculated on the basis of 01(P,, x) and 92(P2, г) 
then value 0:(Р,, х) corresponding to P; (P, < Р; < P, + 0.1) and found according 
to the above approximation are charged with an error smaller than 0.2596. 

The above degrees of the accuracy of the 0;(Р;, zi) determination are satisfactory 
from the practical point of view. 


5. CASE OF A CYLINDRICAL ELECTRODE 


Let us consider an element of an infinitely long cylindrical electrode having 
radius a and length J, see Figure 3a. The solution of appriopriate equation, corre- 
sponding to equation (3), gives 


AJ. n r 


nid s (13) 


F(e,r) - Е (є, a+ d) = 
where AJ, is the diffusive current of electron group Р(=Е)ДЕ. On the basis of formula 
(13), the relationship К! (=) = fi[Fo(s)] cannot be simply determined under the 
condition that F(s,r = oo) = Fo(e), because F(e) — oo аз г — оо. It is the 
principal difficulty of the use of Swift theory for cylindrical electrodes. 

In report [10], the above problem was solved in a rough manner, where it was as- 
sumed that the Swift perturbation is located at the finite region a+d+ several Ae. 
However, this assumption requires grounds and it also requires an additional in- 
formation about the plasma-electrode state. Moreover this procedure can be used 
only in case of the slim electrode, (a << l, А << l). This problem can be also so- 
lved by introduction of the , spherization" of the cylindrical electrode having finite 


length. Let us consider the cylindrical electrode having radius a and length l, see 
Figure 3b. Let us assume that the barrel-shaped surface, placed at distance z from 
the electrode surface and having the area S, = 27|а2 + (a + т)1+ 2z? + raz] is the 
geometric locus of points having identical plasma parameters. Surface S, tends to 
the sphere as r tends to the infinity. Thus we can write 


d|F(s)Ae| | 
dz 


AJ, = e2n[a? + (a + z)l+ 22? + naz} (14) 


Resolving equation (14) in limits d+ «, we obtain 


AJ, 1 1 i, 44+ ra Le VA 


кт deme var 


where A = z(a + pec 8(a? — al). 

Putting Р, (є) in formula (6) and integrating over interval У + oo, we obtain 
the expression for J,(V) which is identical to expression (8). The new formula for 
J.(V) has the new coefficients Ye, Р., determined аз: 


>< n5. За? t al 4d + za +1+УА 
2+4: > 2 мул 4d+ra+l- VA 
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The further use of the Swift theory for the cylindrical probe is the same as for the 
spherical probe. 


6. EXAMPLE OF THE USE OF THE EXTENDED SWIFT THEORY 


Herrmann et al. [11] presented the results of the probe measurement in the 
neon plasma for which the Druyvesteinian distribution is particularly characteristic. 
The electron concentration determined with the cylindrical probe by Herrmann was 
twice bigher than that determined with the spherical probe. This discrepancy could 
be caused by small size of the spherical probe which was 0.6 mm in diameter. The 
slim cylindrical probe has the good contact with the plasma. The small spherical 
probe practically collects the electron current only with half of its surface. Another 
half is shaded by the probe leg. For the better result one should use а probe of a 
larger size. | 

Under the conditions of the Herrmann experiment (P = 1.29 4orr) the mean 
free path of the electrons was 2:1 mm. If the probe of 1.2 mm in diameter is used, 


the ratio i equals 0.6. This value of — is the limit value permissible by Swift in 


his theory. 1a this case and at higher pressures the correcting function, determined 
with integral 0(Р, z), ought to be used. Figure 4 gives the example of the use of 
the-extended Swift theory. 

It is seen that the discrepancy of the characteristics determined with the funda- 
mental theory and those determined with the extended Swift theory is considerable 


if - is high. Moreover it can simply be deduced that the density of the electrode 
А 


е 
electron current decreases when the electrode radius increases. At higher gas pres- 
sures (several torr) the use of the extended Swift theory is necessary even for the 
small electrode as that having 1.2 mm in diameter. 


7T. CONCLUSION 


The generalization of the Swift theory permits to use it at every values of the 
ratio of a electron mean free path and a probe radius. This generalization can be 
used for typical electron energy distributions: Maxwellian and Druyvesteinian. One 
should bé stressed that the Swift theory is important because it takes to account 
the local impoverishment of the plasma occurring around the electrode. 

Swift [3] has given one’s attention to determine of the electrode potential in 
relation to the plasma. The results of this work can simply be incorporated in the 
Swift’s considerations relative to the potential determination. 

In order to correct the energy distributions, measured with an electrostatic 
probe, we have defined a new special integral. This integral can be useful in the 
solving of the other problems. In regards to the physical character of this report 
only the fundamental properties of the above integral are shown. 
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Fig. 2. Integral 6(P, z) as function of parameters P and arguments 2 
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Fig. 4 The current-voltage characteristics of the spherical probe placed in the Druyvesteinian 
plasma. The characteristics were calculated according to the extended Swift theory for several 
values of ratio а/Ле: a — 2$ 0, b — 1, c — 2, d,— 5. Electron current J(V) was normalized 


A gSnav T 
ELO 


According to Swift [3], d = Àe was assumed 


to value J (0) forecasted on the basis of the fundamental Druyvestein theory. 
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Table 1. Values of integrals 6(P, z)* Б 


0.3 0.4 
: 2.1073 0 | 2.2075 0 | 2.3229 0: 2.4581 0 
0.2 | 1.4157 0 | 1.4976 01 1.5952 0 | 1.7111 0 
0.3 | 1.0190 О į 1.0881 0/| 1.1694 0! 1.2669 0 
0.4 | 7.5267 -1 | 8.0875 -1 | 8.7527 -1 | 9.5584 -1 
0.5 | 5.6033 -1 | 6.0518 -1 | 6.5877 -i | 7.2422 -1 
0.6 | 4.1714 4.5249 4.9497 -1 | 5.4725 -1 
0.7 | 3.0913 -1 | 3.3655 -1 | 3.6972 -1 | 4.1079 -1 
0.8 | 2.2739 -1 | 2.4834 -1 2.7382 -1 | 3.0555 -1 
0.9 | 1.6568 -1 | 1.8145 -1 | 2.0070 -1 | 2.2481 -i 
1.0 | 1.1941 1.3108 .-1 | 1.4540 -1 | 1.6342 -1 
1.1 | 8.5023 9.3536 -2 | 1.0401 -1 | 1.1725 -1 
1.2 | 5.9756 -2 | 6.5863 -2 | 7.3403 -2 8.2006 -2 
1.3 | 4.1425 -2 | 4.5734 -2 | 5.1071 -2 | 5.7863 -2 
1.4 | 2.8307 -2 | 3.1299 -2 | 3.5013 -2' 3.9756 -2 


1.9058 


2.1100 «2 | 2.3642 -2 | 2.6897 -2 


1.6 | 1.2636 1.4007 -2 | 1.5717 -2 | 1.7912 -2 
1.7 | 8.2479 -3 | 9.1523 -3 | 1.0282 -2 | 1.1737 -2 
1.8 | 5.2978 -3 | 5.8844 -3 | 6.6190 -3 | 7.5667 -3 
1.9 | 3.3477 -3 | 3.7217 -3 | 4.1908 -3 | 4.7971 -3 
2.0 | 2.0806 -3 | 2.3149 -3 | 2.6092 -3 | 2.9903 -3 
2.1 | 1.2715 -2 14157 -3 | 1.5971 -3 | 1.8324 -3 

1.1036 -3 
2.3 | 7.6394 -4 | 5.0291 -4 | 5.6823 -4 | 6.5322 -4 
2.4 | 2.6180 -4 | 2.9201 -4 | 3.3016 -4 | 3.7987 -4 


2.5 | 1.4928 -4 | 1.6659 -4 | 1.8848 -4 | 2.1702 | -4 
9.3370 -5 | 1.0569 -4 | 1.2179 -4 
2.7 | 4.6019 -5 | 5.1402 -5 | 5.8218 -5 | 6.7130 -5 
2.8 | 2.4872 -5 | 2.7793 -5 | 3.1494 -5 -5 
2.9 | 1.3202 -5 | 1.4757 -5 | 1.6730 -5 | 1.9314 -5 
3.0 | 6.8812 -6 | 7.6947 -6 | 8.7269 -6 | 1.0080 -5 


* The reading manner of the 0(P, х) values is given in the following examples: 


2.2 | 7.6394 JE: -4 | 9.6096 -4 
| 
2.6 | 8.3628 55 | 

| 


d 
e 
© 
e 
= 


6(0.1,0.1) = 2.1073 -10° — 6(0.9,3.0) = 4.6864 . 1075. 
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0.8 
3.4688 
2.6157 
2.0581 
1.6334 
1.2930 


8.2508 -1 


| 
6.1325 -1 | 7.0108 -1| | 1.0154 
4.6324 -1 | 5.3316 -1 | 6.3242 -1 | 7.8898 
3.4635 -1 | 4.0117 -1 | 4.7978 -1 | 6.0539 
2.5601 -1 | 2.9823 -1 | 3.5933 -1 | 4.5807 
1.8668 -1 | 2.1880 -1 | 2.6541 -1 | 3.4158 
| 
1.3456 -1 | 1.5829 -1 1.9320 -1 | 2.5084 
9.5530 -2 | 1.1285 -1 | 1.3853 -1 | 1.8132 
6.6825 -2 | 7.9250 -2 | 9.7779 -2 | 1.2895 
4.6038 -2 | 5.4788 -2 | 6.7920 -2 | 9.0208 
3.1223 -2 | 3.7277 -2 | 4.6415 -2 | 6.2048 
2.0839 -2 | 2.4953 -2 | 3.1195 -2 | 4.1957 
1.3684 -2 | 1.6429 -2 | 2.0615 -2 | 2.7884 
8.8377 -3 | 1.0636 -2 | 1.3393 -2 | 1.8211 
5.6124 -3 | 6.7701 -3 | 8.5519 -3 | 1.1686 
3.5040 -3 | 4.2355 -3 | 5.3662 -3 | 7.3668 
| 

2.1503 -3 | 2.6042 -3 | 3.8084 -3 | 4.5615 
1.2968 -3 | 1.5733 -3 | 2.0039 -3 | 2.7741 
7.6846 -4 | 9.3388 -4 | 11922 -3 | 1.6568 
4.4739 -4 | 5.4455 -4 | 6.9672 -4 | 9.7162 
2.5586 . -4 | 3.1185 -4 | 3.9984 -4 | 5.5947 
1.4372 -4 | 175408 -4 | 2.2533 -4 | 3.1627 
7.9289 -5 | 9.6884 -5 | 9.2468 -5 | 1.7552 
4.2954 -5 | 5.2546 -5 | 6.7737 -5 9.5623 
2.2849 -5 | 2.7980 -5 | 3.6133 -5 | 5.1133 
1.1933 -5 | 14628 -5 


1.8914 -5 | 2.6837 


0.9 
4.1329 
3.2345 
2.6206 
2.1331 
1.7272 


1.3848 
1.0966 
8.5664 
6.5688 
4.9899 


- 3.7180 


2.7244 
1.9627 
1.3896 


9.6688 


6.6089 
4.4373 


| 2.9261 


1.8949 
1.2048 


7.5221 
4.6102 
2.7737 
1.6380 
9.4947 


5.4014 
3.0156 
1.6522 
8.8831 
4.6864 


(2) ©) ie (=) 
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8. APPENDIX 
exp —t? 
(Е -— T TS 


f* exp —t? _ exp—-z? -[ ume 
J, (t-Pzp^ 2201- P jJ, @- Pz)? 


(i) The integral Г dt of the expression (9) was calculated as follows: 


а EXD 2 С Я f? exp —t? 
_ 222(1- Р)? J, (t Pp TAE (t — Pz)? 


The next use of the integration by parts and of the substitution t — Pz — z gives 


f° ехр 12, 1 p : 3a 
J, @- Psp = lin» pep] (ex s 


+2Pz erfc (т) — (1— 2P?z?)86(P, z) 
(ii) According to the above scheme, the integral defined as 


20 t? exp —t? 
OP, =) въ - | (Е рар 


can be transformed to the expression containing integrals егіс (т) and 6(P, т). 
р, v are natural numbers and y < v. 
(iii) Putting t? = z, the integral 0(0, х) can be evaluated as 


коа) = [^ set = dt = 5822). 


Jr 
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STRESZCZENIE 
! 

W świetle teorii Druyvesteina i Swifta badano prąd elektronowy sondy. Wprowadzenie nowej 
funkcji specjainej pozwala okreslić prąd elektronowy sondy dla dowolnego stosunku promienia 
sondy do średniej drogi swobodnej elektronów. Wartości tej funkcji, obliczone numerycznie sa 
zestawione w tabeli. Przedstawiono zastosowanie uogólnionej teorii Swifta dla sondy sferycznej. 
Ponadto pokazano, ze nowa funkcja specjalna jest elementem rozwiązania innych nieelementarnych 


calek. 


PE3IO ME 
B рамках теории Дрювестейна и Свифта исследовался электронный ток зонда. 
Введение новой функдии позволяет определить электронный ток зонда при лю- 
бом отношении радиуса зонда к свободному пробегу электронов. Значения 
этой функции, вычисленные с помощью ЭВМ, указаны в таблице. Показано 
применение обобщённой теории Свифта к сферическому зонду. Показано также 


использование новой специальной функции для решения других невлементарных 
интегралов. 


Zloione 31.У 1.1987 


